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The formation and evolution of galaxies is imprinted on their stellar population radial gradients. Two recent articles present 
conflicting results concerning the mass dependence of the metallicity gradients for early-type dwarf galaxies. On one side, 
Spolaor et al. show a tight positive correlation between the total metallicity [Z/H] and the mass. On the other side, in 
a distinct sample**, we do not find any trend involving [Fe/H] (Koleva et al.). In order to investigate the origin of the 
discrepancy, we examine various factors that may affect the determination of the gradients: namely the sky subtraction and 
the signal-to-noise ratio. We conclude that our detection of gradients are well above the possible analysis biases. Then, 
we measured the [Mg/Fe] relative abundance profile and found moderate gradients. The derived [Z/H] gradients scatter 
around -0.4dex/r e . The two samples contain the same types of objects and the reason of the disagreement is still not 
understood. 
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1 Introduction 

The stellar populations of galaxies hold a fossil record of 
their formation and evolution. Their study may allow to 
discriminate between the various formation scenarios pre- 
dicting different spatial gradients and different relations be- 
tween the gradients and the mass of the galaxy. 

In the classical monolithic scenario, galaxies formed in 
a dissipative collapse (Larson, 1974; Arimoto & Yoshii, 
1987), where stars remain on their orbits and do not mix. 
The gas, enriched in metals from the evolved stars, flows 
into the centre, generating negative metallicity gradients 
(i.e. higher metal content in the centre than in the outskirts). 
The chemo-dynamical simulations Matteucci & Tornambe 
(1987, see also Kawata & Gibson, 2003) predict, that due 
to their small potential, dwarf galaxies, are unable to retain 
metals and thus their gradients are shallower or almost inex- 
istent. Therefore relations between the mass and the metal- 
licity or the metallicity gradient are expected. 

A competing (or complementary) formation scenario is 
the hierarchical clustering (e.g. Cole et al., 1994), where 
small dark matter halos merge to produce bigger galaxies. 
It may be thought that these violent events mix the stellar 
population and erase the gradients (e.g. White, 1980), but 
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** Based on observations made with ESO telescopes at La Silla Paranal 
observatory under program ID076.B-0196 



some simulations (e.g. di Matteo et al., 2009), at variance, 
show a gradient preservation, due to the violent relaxation 
which does not mix the orbits. 

The interpretation of the metallicity gradients is not 
straightforward. Nevertheless, gradients observation could 
constrain model's predictions. Several studies on big sam- 
ples of giant elliptical galaxies have been already performed 
(e.g. Sanchez-Blazquez et al., 2006). However, due to their 
low-surface brightness, the dwarf galaxies, are still very 
much 'terra incognita'. 

Recently we published a paper on the stellar content 
of 16 dwarf galaxies in the Fornax cluster and in groups 
(Koleva et al., 2009a). One of our main conclusions was 
that most of our galaxies hold a strong negative metallicity 
gradient. The star formation history analyses revealed that 
these gradients are built at early epochs. We did not find 
any correlation between our gradients and the velocity dis- 
persion (proxy of the mass) in our sample (Fig.l). However, 
this result is challenged by a recent letter from Spolaor et al. 
(2009), using a sample of 14 dwarfs in the Virgo and For- 
nax clusters observed with GMOS at the Gemini South tele- 
scope in medium resolution, long-slit, optical spectroscopy. 

To analyse their data, Spolaor et al. (2009) binned 
them along the slit direction to achieve a minimum signal- 
to-noise of 30 @ 5200 A. They measured Lick indices 
(Worthey et al., 1994; Worthey & Ottaviani, 1997) and de- 
rived ages and chemical compositions by comparing to the 
Thomas et al. (2004) single stellar population (SSP) models, 
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using their own method (Proctor & Sansom, 2002). They 
estimated the gradients performing a linear least-squares 
fit to the radial metallicity profile weighted by their er- 
rors. The fit was done in logarithmic space of the radius, 
A[Z/H]/A\og(r/r e ), where r varies from 1 arcsec to one 
effective radius (r e ). 



In Sect.2 we present the sample, observations and anal- 
ysis. In Sect.3.1 we investigate the systematics (signal-to- 
noise ratio, bad sky-subtraction) which may affect the de- 
rived metallicity gradients. As Spolaor et al. used the total 
metallicity [Z/H]we also measured [Mg/Fe] and discuss 
the differences between A[Fe/H] and A[Z/H] in Sect.3.2. 
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Fig. 1 Metallicity gradients vs. a and M# (top and bottom 
respectively) of the 2 samples of dwarf elliptical galaxies, 
one from Spolaor et al. (2009, blue, open squares) and other 
is this work (filled symbols). Green circles are galaxies from 
Vazdekis et al. (2001), red circles from Koleva et al. (2009a) 
and a red star for NGC 205(Simien & Prugniel, 2002). The 
names of our galaxies are also plotted. 

There is no galaxy in common between the two samples 
and the causes of the disagreement may be multiple: choice 
of the sample, data reduction or analysis issue. The primary 
purpose of the present paper is to explore these possibili- 
ties. In order to directly compare our results with Spolaor 
et al. (2009) we will determine our metallicity gradients in 
an identical way. 



2 Data and analysis 

For the purpose of the present paper we collected 20 early- 
type galaxies (Table 1) with velocity dispersions less than 
150 km s -1 from three different data sets (VLT, OHP and 
WHT). We analysed them with the same method. 

2.1 OHP data 

We used long slit data of NGC 205 (Simien & Prugniel, 
2002), which were taken with CARELEC spectrograph on 
1 .93 m telescope in Observatoire de Haute-Provence (OHP). 
The observations were carried from January 2001 to Jan- 
uary 2003 and are partly used in Simien & Prugniel (2002) 
to derive the internal kinematics. The slit width was set 
to 1.5", which results in an instrumental broadening of 
R = 5050 ( FWHM, or a ms w 25 km s" 1 ). The wavelength 
range covers AA 4700 - 5600 A. NGC 205 is an emblematic 
galaxy for the dwarf elliptical class. 

2.2 VLT data 

We obtained long slit spectra of 16 dwarfs elliptical galax- 
ies (dE) with FORSl,2 mounted on the VLT. These galax- 
ies are situated in the Fornax cluster and the NGC 5044, 
NGC 5898 and NGC 3258 groups, i.e. in different environ- 
ments. This sample covers a diversity of dwarfs: nucleated, 
non-nucleated, with or without spiral substructures, with or 
without gas (Koleva et al., 2009a, tablel). These dwarfs 
have Mb between -17.8 and -15.4, while their velocity dis- 
persion 1 ranges from 39.5 to 87.0kms~ 1 . In other words 
they are similar to the dEs prototype - NGC 205. 

The group galaxies were observed with FORS2 with 
the GRIS1200g+96 grism, resulting in AA 3300-6200 A and 
fms = 74kms _1 . The Fornax data were observed with the 
GRIS 600B+22 grism on FORS 1, which results in AA 4335- 
5640 A and a ins = 64 km s" 1 at 5200 A. 

2.3 VAKU sample 

We also used three of the 6 Virgo cluster galaxies analysed 
in Vazdekis et al. (2001). These three objects are situated 
in the transition region between normal and dwarf ellipti- 
cal galaxies (rr » 120km s~ x and Mb ~ -18). The long slit 
observations were made with the William Herschel Tele- 
scope (WHT) in La Palma. The wavelength coverage is 
AA « 4000 — 5000 A, and the instrumental resolution is 

2.4 A, or R« 2000. 



1 luminosity weighted mean velocity dispersion from De Rijcke et al. 
2005 
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Table 1 Properties of our galaxies sample. In column 1 we list the galaxies' names, in columns 2 and 3 the equatorial 
coordinates, in column 4 the B-band apparent magnitude trig, in column 5 the Schlegel galactic extinction (Schlegel et al., 
1998, ABG,PLElNPOT), in column 6 the radial velocities (measured with ULySS), in column 7 the k-correction (from 
COCOR, PLEINPOT, using the listed cz values), in column 8 the absolute magnitudes Mb in B-band (computed as described 
in the text), in column 9 the effective radius derived from the /?-band images (De Rijcke et al., 2005), in column 10 
the luminosity-weighted velocity dispersion (De Rijcke et al., 2005), in column 11 the metallicity gradients (derived as 
described in the text). 
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a from HyperLeda database (Prugniel & Simien, 1996); 6 from Simien & Prugniel (2002); c from Koleva (2009); d from Vazdekis et al. 
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2.4 Analysis 

The SSP-equivalent profiles (see Fig. 1, 2 of Koleva et al. 
2009a) were obtained by binning the 2D spectrum in the 
spatial direction to achieve a signal-to-noise ratio of 20 and 
by analysing the resulting ID spectra with ULySS 2 (Univer- 
site de Lyon Spectroscopic analysis Software, Koleva et al., 
2009b). 

ULySS is designed to fit any linear combination of non- 
linear components (weighted by W) against a spectrum 
fobs (A). A model constructed as such, can be convolved 
with a line-of-sight velocity distribution (LOS VD) and mul- 
tiplied by a n th order Legendre polynomial, P n (X).: 

F obs (X) - P n (A) x \ LOSVD(v sys ,a,h3,M) 

i—k \ 

®J2 W * CMP 1 (oi,a 2 ,...,A)]. 

(1) 

The LOS VD is a function of the systemic velocity, v sys and 
the velocity dispersion a and may include a Gauss-Hermit 
expansion (/i3 and M, van der Marel & Franx, 1993). A 
is the logarithm of the wavelength (the logarithmic scale is 

2 http://ulyss.univ-lyonl.fr 



required to express the effect of the LOSVD as a convolu- 
tion). In this particular case, the CMP (i = 1) is a SSP 
constructed with Pegase.HR stellar population models (Le 
Borgne et al., 2004) using Elodie.3.1 empirical stellar li- 
brary (Prugniel & Soubiran, 2001; Prugniel et al., 2007b). 
Consequently, the parameters a\ , 0,2 are the age and metal- 
licity of the SSP. For the a-elements analysis we added 
another parameter 03 which expresses the [Mg / Fe] abun- 
dance. In the latter case the stellar population models were 
build with a semi-empirical library, produced by combina- 
tion of Elodie.3. 1 and Coelho et al. (2005) libraries, includ- 
ing the non-solar abundance of [Mg/Fe] (Prugniel et al., 
2007a; Koleva et al., 2008). More precisely, Elodie.3.1 
was differentally corrected for [Mg/Fe] = 0.0 and 0.4, us- 
ing Coelho et al. library. The analysis selects the optimal 
[Mg/Fe] value interpolated in this range. 

Some of the galaxies formation models predict relation 
between the mass of a galaxy and its metallicity gradient. 
As a proxy to the mass we used the velocity dispersion 
and the absolute magnitude. We took the velocity disper- 
sion values from the literature (see Table 1). The absolute 
magnitudes were computed using cz, mj (from Tablel) and 
Hubble constant H =70 kms -1 Mpc^ 1 . We corrected for 
Schlegel et al. (1998) galactic extinction and reddening (de- 
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rived with ABG and COCOR respectively from the PLEIN- 
POT 3 package). 

We did not find any obvious correlation between our 
gradients and the mass of the galaxies. From Fig. 1 one can 
notice that if we had data only from Vazdekis et al. (2001) 
and several dwarfs from Koleva et al. (FCC 288, DW1, 
DW2, FS 131) then we could have found a mass-gradient 
relation. We note also that, some of our flattest galax- 
ies (e > 0.3; DW1, DW2, FCC 046, FCC 288, FS029, 
FS 131) follow a trend with M B and a. However, FCC 204, 
which is also fiat (e = 0.61), is an outlier. We cannot ex- 
clude the possibility that the contrasting results are due to 
the small size of the sample or to the actual selection critre- 
ria. 



3 Validation 

3.1 Noise and sky subtraction 

The two main features of the data reduction and analysis 
which can affect the derived metallicity gradients are a pos- 
sible bad sky subtraction and the decreasing signal-to-noise 
from the centre to the outskirts. Our measurements are al- 
most insensitive to flux calibration problems (as the Lick 
indices). The shape of the flux is modeled by a multiplica- 
tive polynomial, during the fitting and cannot be a source of 
biases. 

To investigate how the signal-to-noise and imperfections 
of the sky subtraction can affect the measured metallici- 
ties, we produced a ID model of the centre of FS 373, us- 
ing the Pegase.HR population synthesis code. The choice 
of this particular galaxy, was driven by the its strong metal- 
licity gradient (-0.46 dex). To construct the model, we use 
the SSP-equivalent parameters as derived from our analysis, 
i.e. age= 1.55 Gyr and [Fe/H] = 0.0 dex. We broadened the 
model to the velocity dispersion (instrumental plus physi- 
cal) of the observed spectrum. 

To investigate the metallicity dependence on the S/N we 
added Gaussian noise to the model using SNR and NSIMUL 
keywords in the call of ULySS. The results are shown in 
Fig. 2. The errors become larger when the S/N decreases, 
an expected behaviour. However, we do not observe any im- 
portant systematics in the measured metallicity values down 
to S/N =10. 

The sky subtraction, especially for diffuse objects like 
the dwarf galaxies, is a particularly critical point. To ob- 
tain a negative metallicity gradient, one needs to system- 
atically under-subtract the background. To explore this ef- 
fect, we took an UVES sky spectrum (Hanuschik , 2003) 
and broadened it to the instrumental resolution of FORS2, 
74kms~ 1 , and scaled it to the typical brightness of the 
dark sky. Taking into account the redshift of the galaxy we 
added the two spectra (sky plus galaxy). We simulated an 
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Fig. 2 Dependence of the ULySS metallicity determina- 
tion on the signal-to-noise of the data. With grey line we plot 
the input metallicity and with blue filled circles the output 
values of [Fe /H] . 



under- subtraction of the sky from fully subtracted to un- 
subtracted. The results are shown in Fig. 3. To decrease the 
gradients by 0.2 dex, we need to under-subtract 50 % of the 
sky. For FS 373, where the gradient is -0.46 dex, we would 
need to under-estimate the sky by more than 100%! We 
could possibly make an error on the sky subtraction of 10 % 
in improbable cases, but certainly not as a systematic under- 
subtraction. Our strong gradients cannot be an artifact of a 
bad sky subtraction. At least not in the way we modeled it 
here. 
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Fig. 3 Dependence of the metallicity determination with 
ULySS on the under-subtracted sky in fraction. The grey 
line marks the input metallicity of the model, while the blue 
dots are the derived metallicity. In abscisse, fraction = 1 cor- 
responds to the non subtraction of the sky, and fraction = 
to the exact subtraction. 
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3.2 Iron against metallicity abundance 



4 Discussions 



The metallicity of the stellar population models are labeled 
in [Fe/H] (because the stellar libraries are indexed in this 
manner). Therefore, although the stellar evolution depends 
on the detailed abundance (Salasnich et al., 2000), the full- 
spectrum fitting measures [Fe/H] rather than [Z/H] (the 
differences in the stellar evolution mostly translates as bi- 
ases on the age). 

Spolaor et al. using Thomas et al. 2004 models, derived 
their [Z/H] as a sum of Fe-like plus Mg-like elements (C, 
N, O, Mg, Na, Si). Consequently, if the populations are not 
scaled-solar, the results may not be consistent with ours. 

To make the comparison rigorous, we also derived 
[Z/H] by measuring [Mg/Fe]. We explore the possibility 
that the [Mg/Fe] gradients compensate [Fe/H] to result 
in flat [Z/H] profiles. To test the difference in the mea- 
sured quantities we used our a-enhanced models (Prug- 
niel et al., 2007a; Koleva et al., 2008) as explained in 
Sect. 2. We found in general negative [Mg/Fe] gradients 
(A[Mg/Fe] = -0.10 ± 0.09dex ). Only four galax- 
ies (DW1, FS76, FS 131, FCC 043) have flat or positive 
[Mg/Fe] gradient. We computed the total metallicity using 
[Z/H] = [Fe/H]+0.9S [Mg/Fe] (Thomas et al., 2004). 
We find that the total metallicity gradients (see Fig.4) are 
stronger than the [Fe/H] gradients, and the difference in 
the metallicity definitions cannot explain the difference be- 
tween the two data sets. 
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We investigated the discrepancy in the metallicity gradients- 
mass relation between two recent works, namely Spolaor 
et al. (2009) and Koleva et al. (2009a). While the former 
work shows a very strong and tight relation between the to- 
tal metallicity gradient in low-mass galaxies and their cen- 
tral velocity dispersion or Mb, the latter do not find a hint of 
such a relation. 

As first guess for this discrepancy we pointed the small 
size of the samples and different selection criteria. Such a 
scenario is possible, since some of our galaxies fall on Spo- 
laor et al. relation. This hypothesis may be supported by 
the fact that the Fornax galaxies of Spolaor et al. (2009) 
are in general flat (e > 0.3), and we have shown that the 
gradients in such galaxies are usually shallower. However, 
in their Virgo sample all but one galaxy have e < 0.3 4 . 
Hence, this explanation is ruled out and flattening cannot be 
a reason for our disagreement. 

Another possible explanation would have been that 
while the dwarf galaxies exhibit a strong negative 
A[Fe/H], they possess the same [a/Fe] gradient but with 
a positive sign. Thus, measuring the total metallicity as a 
combination of [Fe/H] and [a/Fe], we would have a total 
A[Z/ H]=Q. We measured the [Mg/Fe] gradients in our 
sample using unpublished semi-empirical models (Prugniel 
et al., 2007a). These models are preliminary but the relative 
measurements of [Mg/Fe] are certainly reliable. Opposite 
to the expectations we found that (most of) our galaxies pos- 
sess negative [Mg/Fe] gradients. Consequently, the total 
metallicity gradients are stronger than the [Fe/H] gradi- 
ents. 

Negative iron and [Mg/Fe] gradients will imply that, 
in the early epochs, when the feedback is a strong regu- 
lator, the star formation starts in the centre, where the gas 
is the densest. After the gas is exhausted, the star forma- 
tion stops first in the centre and continues longer at larger 
radius. This process forming a negative [Mg/Fe] gradient. 
To form the negative [Fe /H] gradient we need that cool, en- 
riched gas flows again into the centre, triggering star forma- 
tion. In other words, having repeating episodes of star for- 
mation. This scenario is in agreement with recent SPH sim- 
ulations on the formation of dwarf galaxies (Valcke et al., 
2008, fig. 10) 

We also investigated the effect of a bad sky-subtraction 
and low signal-to-noise on our metallicity measurements. 
We showed that our metallicity measurements are not biased 
down to S/N =10 and that to have zero metallicity gradients 
we need to under-subtract the sky with more than 50 %. We 
concluded that these two effects cannot be the cause of our 
strong gradients. 

As a conclusion, according to the tests performed in this 
paper, the origin of the differences between the two articles 
remains unknown. 



We used the ellipticity measurements from the HyperLeda database 
(Paturel et al., 2003), http : //leda . univ-lyonl . f r 
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